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We examine the dynamics of a qubit stored in the motional degrees of freedom of an ultra-cold 
ion in an ion trap which is subject to the decoherence effects of a finite-temperature bath. We 
discover an encoding of the qubit, in two of the motional modes of the ion, which is stable against 
the occurrence of either none or one quantum jump. For the case of a zero-temperature bath we 
describe how to transfer only the information concerning the occurence of quantum jumps and their 
types to a measuring apparatus, without affecting the ion's motional state significantly. We then 
describe how to generate a unitary restoration of the qubit given the jump information, through 
Raman processes generated by a series of laser pulses. 
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I. INTRODUCTION 



Is it possible to preserve quantum coherence in systems 
that are irreversibly coupled to the environment? Apart 
from being of fundamental interest, this question is of 
particular significance in the area of quantum informa- 
tion processing jjj , where the first two-bit quantum gate 
operations have been demonstrated in quantum optical 
systems such as trapped ions , and optical cavities j| . 
In this context the above question has been answered 
with "yes" through the developement of quantum error 
correction concepts, in which quantum coherence is pre- 
served by exploiting multi-particle entanglement . The 
realization of quantum gates using nuclear magnetic res- 
onance and spin-spin coupling within a molecule Q , has 
already allowed a test of simple ideas of such correction 
concepts ||. 

At present, perhaps the most promising implementa- 
tion of a quantum information processor which comprises 
a small number of qubits is based on a linear ion trap 
as first proposed by Cirac and Zoller In most cur- 
rent formulations of such an ion trap computer ^,|) , the 
qubits are stored in two electronic ground or metastable 
states of individual ions while the quantum "bus" con- 
necting the various qubits is mediated by the ion's col- 
lective motional state either through the center-of-mass 
(CM) or higher excited modes. At the same time, recent 
experiments have shown that the decoherence effects in 
systems of trapped ions are dominated by decoherence 
of the motional state |^, PJ0|. This is also supported by 
theoretical investigations |11H15|1 . 

The aim of this paper is to explore possibilities for 
motional error correction and to present a scheme which 
may be used to preserve quantum coherences in the 



motional state of trapped ions. We consider a single 
trapped ion whose motional degrees of freedom undergo 
finite-temperature dissipative dy namics as described by 
a thermal master equation In this situation we 

describe how quantum information, when stored in the 
motional (bosonic) degrees of freedom of a trapped ion, 
can be encoded in a way which allows its active sta- 
bilisation. For the zero-temperature case we explicitely 
demonstrate how such a stabilisation can be implemented 
in a trapped-ion system. This stabilisation can correct for 
the loss of a single excitation in the motional state and 
needs only a single ion while making use of two of the 
three degrees of freedom of the ion's motion. Such sta- 
bilisation techniques could be used to stabilise the "bus" 
from the effects of decoherence or, more speculatively, 
the stored qubits in a bosonic quantum computer ]l5[ . 

We will describe decoherence using the quantum jump 
approach |l6| and base our stabilisation scheme on the 
concept of "quantum jump inversion" [ p~7| . To begin 
with, we briefly review some of the underlying ideas and 
give an overview of the method that we will use to sta- 
bilise the quantum information. In Section |l| we out- 
line the physics of the motion of the trapped ion while 
in Section 111 we go into detail regarding the short-time 
dissipative evolution and the derivation of the specific 
qubit encoding which allows the active stabilisation of the 
stored quantum information. In Section IV we describe 



our method of periodically detecting whether a quantum 
jump has occurred and if so, determine the type of jump. 
This step will involve a projective measurement of the 
ion's electronic states which does not significantly dis- 
turb the motion, a process which we discuss in Section 
IV B|. Finally, in Section [v|, we show how the unitary 



inversion of the detected decoherence effects (quantum 
jumps) can be effected through adiabatic processes. 
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A. Decoherence and Quantum Jump Inversion 
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Our scheme is based on the concept of "quantum jump 
inversion" which was first presented by Mabuchi and 
Zoller jl7| and represents a special example of the general 
formalism of reversible quantum operations as studied in 
detail by Nielsen and Caves [Q . We will make some use 
of their results in the following and briefly review them 
here, together with the underlying description of deco- 
herence through the formalism of quantum trajectories 
Jl9| , p0[ and the quantum jump approach [jl6| . 

In quantum optics, decoherence is described by cou- 
pling the system to an infinitely large reservoir which 
models the environment. The couplings between the sys- 
tem operators and the reservoir operators are usually as- 
sumed to be weak and linear. This allows one to derive a 
Markovian master equation for the dynamics of the sys- 
tem reduced density operator p , after tracing over the 
reservoir degrees of freedom Km. The general form of 
this master equation is 

N 

~ E 2 { d ^mP + P& m C m - 2C mP Cl) , (1) 

m— 

where H is the system Hamiltonian and the operators C m 
describe the different decay channels to the environment. 
The C m originate from the linear coupling between sys- 
tem and environment and are proportional to the system 
operators which couple to the reservoir. 

An alternative approach to open quantum systems has 
been developed which describes the system through an 
ensemble of pure state evolutions through the formalism 
of quantum trajectories Jll||2(|. In this method the sys- 
tem is continuously monitored by performing measure- 
ments on the state of the reservoir. Different measure- 
ment strategies lead to different conditioned evolutions of 
the system which are termed unravelings. More specif- 
ically, in the quantum jump method pq| , one obtains 
an unraveling of the master equation as arising from the 
conditioned evolution of the system where the condition- 
ing is supplied by detecting the loss of system excitations 
through the counting of quanta in the environment. In 
the absence of any measurement counts the system evo- 
lution is governed by the non-Hermitian Hamiltonian 

h N 

H eff =H~ l -J2clC™, (2) 

m— 

which differs from the system Hamiltonian H due to 
the information that we acquire from not observing any 
counts. On the other hand, the information obtained 
from a detector count at a time t, associated with the 
decay channel to, conditions the state of the system ac- 
cording to 



where \ip(t— )) and \ip(t+)) describe the state of the sys- 
tem immediately before and after the quantum jump, 
respectively. 

In their article, Mabuchi and Zoller have formu- 
lated conditions for the reversibility of quantum jumps 
through unitary operations. We adopt here the more 
general notation of Nielsen and Caves for the re- 
versibility of quantum operations. Nielsen and Caves 
have shown that, given the initial state \tp) of a quantum 
system is known to lie in a specific subspace Ho C H, of 
the entire system Hilbert space H, and the subspace Ho 
has the property 

= M 2 , V|V) G Ho , (4) 

then the action of the quantum operation as described 
by A is unitarily reversible. In other words, given 
Eq. (|4|) then there exists a unitary operation U, such that 
UA\tp) — p\ip) , V|V) G Wo, where p is a real constant [ pj| . 
To illustrate the condition in Eq. (|j), consider a system 
with a single bosonic degree of freedom, such as the elec- 
tric field inside a single-mode resonator, and where Ho 
is a two dimensional subspace of H and is spanned by 
two basis states and \ip-). If the quantum operation 
A = a, and describes a quantum jump corresponding to 
the loss of a photon from the resonator, then Eq. (|J) im- 
plies that all states € Ho, possess the same mean 
cavity excitation number. Alternatively, Eq. (^J) is sat- 
isfied if the basis states \ip+) and \ip-) are orthogonal, 
where \ip±) possess the same mean excitation number 
and remain orthogonal after the quantum jump. It is 
clear that in this situation the coherence present in the 
subspace Ho (in the form of a coherent superposition of 
the states \ip+) and |^>_)) is not damaged through the 
action of the quantum operation A and it can be re- 
stored to the state prior to the quantum jump through 
a unitary process. Mabuchi and Zoller have exploited 
this fact to preserve quantum coherences in a specially 
constructed situation in cavity quantum electrodynam- 
ics (CQED) fTJ. Following their work, Mensky |2l) has 
studied the inversion of quantum jumps under the aspect 
of measurement reversibility. More recently, Vitali et al. 
p2~f have proposed CQED schemes which are similar in 
spirit to that of Ref. (l7]] to protect quantum states in 
cavities employing continuous feedback methods. 



B. Stabilization Strategy 

In the following we outline how the strategy of quan- 
tum jump inversion can be modified and implemented to 
stabilise a qubit which is stored in the motional degrees 
of freedom of a trapped ion. The implementation of a 
quantum jump inversion scheme requires a priori knowl- 
edge about the decay channels that describe the effects 
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of decoherence as in 

Eq.(@). We will assume that the 
motional degrees of freedom undergo finite-temperature 
dissipative dynamics which are described by a thermal 
master equation As we will show below, to derive a 
qubit encoding which allows the unitary inversion of the 
quantum jumps associated with this type of decoherence, 
it is sufficient to examine the case of a zero-temperature 
bath. This will considerably simplify our analysis. Irre- 
spective of the particular choice of bath, the exact inver- 
sion of a quantum jump will require: (i) A special choice 
of motional subspacc Ho in which to encode a qubit such 
that the quantum jumps associated with a given type 
of decoherence are unitarily reversible, (ii) a mechanism 
for detecting quantum jumps without much disturbance 
of the motion, and (iii) a procedure which accomplishes 
the unitary inversion of the detected decoherence effect. 
While (i) is difficult, to satisfy (ii) and (iii) represents a 
serious challenge for an operational quantum jump inver- 
sion strategy. 




FIG. 1. Periodically performed projective measurements 
on the system: With a period r we interrogate the dissipa- 
tive system evolution. Step (i) serves to determine whether 
a quantum jump has occurred in the preceding time interval 
r. If a quantum jump is detected in step (i) this is reversed 
through a unitary operation in step (ii) and the quantum 
coherence in the system is restored to the original state. 

Before we go into more detail about how the above 
mentioned requirements can be met in the physical sit- 
uation of a trapped ion we mark the following general 
considerations. The effects of motional decoherence that 
have been observed in recent trapped ion experiments are 
believed to be due to technical noise fiofl . In this situ- 
ation, even though one can obtain master equations for 
the motional dynamics of the ions I0J12J, it is difficult to 
isolate the particular environment to which the system is 
coupled. Even in those cases where a physical quantum 
bath can be identified, it is extremely difficult to monitor 
the system through the quantum information transferred 
to the reservoir. The energies of the motional excitations 
are too small to be easily detected. Thus, to determine if 
a quantum jump has occurred we propose to periodically 
perform projective measurements on the system directly 
without gaining any information on the stored quantum 
superposition, as illustrated in Fig.^j. This is a signifi- 
cant difference to the original formulation of a quantum 
jump inversion scheme p!7y , where the dissipative system 
dynamics are continuously monitored by imposing detec- 



tors directly on the environment output channels and it 
makes our scheme more akin to quantum error correc- 
tion schemes H|. The main consequence of this modifi- 
cation lies in the fact that our quantum jump detection 
scheme cannot reveal when in the time interval between 
successive projective measurements the quantum jump 
occurred. We will however still be able to unitarily in- 
vert the decoherence effects caused by such a jump. 

To deduce the proper encoding for the qubit we exam- 
ine, in Section |n], the effects of motional decoherence as- 
sociated with a zero-temperature master equation given 
that the initial motional state of the system \ip V ib) lies 
in the subspace Ho C H, which is spanned by two basis 
states \ip±), and where H is the entire motional Hilbert 
space. As mentioned above this is sufficient to obtain 
an encoding which allows the unitary inversion of the 
quantum jumps associated with the more general, finite- 
temperature reservoir. Based on the condition given in 
Eq. (|J) we show that for the decoherence effects to be uni- 
tarily reversible both in the case of a detected quantum 
jump, and in the absence of a detection event, requires 
that the basis states \ip±) comprise of at least two bosonic 
modes. We therefore propose to encode the qubit as 



\tpvib) = C+\lp+) +C-\lp-) 



(5) 



in two of the ion's motional degrees of freedom. In this 
situation, the effects of decoherence cause the system to 
evolve into a mixed state 



P ~PoW> )(^ \ +P2|^W 



(6) 



on a time scale which is short compared with the char- 
acteristic decay time. Here, \ip°) £ Ho, IV^ 1 ) € I'll, and 
|V-' 2 ) € 7^2- Since the projective measurements which ac- 
complish the detection of quantum jumps need to dis- 
criminate between the three parts in the above mixture, 
we require the three subspaces Ho, Hi and H2 to be mu- 
tually orthogonal. These constraints will lead us to the 
choice 

\i>±) = (|4,0) + |0,4)± V2|2,2))/2, 

for the basis motional states where \n x ,n y ) denotes a 
two-dimensional motional Fock state. 

The detection of quantum jumps in the motional state 
of the trapped ion represents step (i) in the schematic de- 
piction of our stabilisation scheme given in Fig.|. This 
step consists of two parts. First, a Raman induced 
process serves to entangle the mixed state given above 
in Eq. @ , with orthogonal electronic states of the ion, 
thereby transferring information about the occurrence of 
quantum jumps into the electronic degrees of freedom. 
Following this e ntanglement operation, which we discuss 
in Section IV A, we perform a measurement of the ion's 
electronic degrees of freedom which projects the state 
from the mixture in Eq. (^) into a pure state. The pro- 
jection coul d, in p rinciple, be done using quantum jump 
techniques |23,E4]. However, in the situation where we 



3 



aim to protect the motional state of an ion from the ef- 
fects of decoherence we cannot tolerate the recoil from 
scattering a large number of photons as this would lead to 
motional excitation and irreversibly change the motional 
state p5| . Instead, for the second st ep of the quantum 
jump detection, we propose in Section IV E| an alternative 
projective measurement which is based on the "photon 
gun" scheme discussed by Gheri et al. (2(| , where the ion 
is surrounded by a low-Q optical cavity. Dependent on 
the ion's electronic state a laser pulse triggers the trans- 
mission of a single-photon state out of the cavity where 
the projection can be accomplished through photodetec- 
tion. In the Lamb-Dickc limit |27],|28| we can practically 
eliminate any motional recoil suffered by the ion while 
still obtaining a signal whether a quantum jump has oc- 
curred or not. In effect, we are coupling in another out- 
put channel to interrogate the system in addition to those 
channels coupling the system to the reservoir. 

The unitary inversion of the detected quantum jumps 
which is schematically depicted as step (ii) in Fig. [I] repre- 
sents perhaps the most essential part of the stabilization 
scheme. This is addressed in Section In the case where 
no quantum jump has been detected the ion is left in its 
initial motional state as given in Eq. ([s]) and no further 
manipulation is required to restore the quantum infor- 
mation to the original basis states. In the case where a 
quantum jump has been detected, we describe, using two 
adiabatic transfer processes and an intermediate resonant 
step, how one can unitarily invert the decoherence effects 
associated with this quantum jump. The involved pro- 
cesses have the advantage that they are mostly achieved 
via adiabatic passage and thus are not subject to noise 
due to spontaneous emission nor do they suffer from the 
strict timing constraints imposed by resonant processes. 
The laser arrangements needed to effect this unitary in- 
version of the detected decoherence effects are not overly 
complex and may be experimentally feasible. 

In our analysis of the processes that constitute step (i) 
and (ii) in the schematic depiction of our stabilisation 
scheme in Fig. [l] we will neglect all effects of decoherence 
during the corresponding operations. In this situation, if 
the detection and inversion processes are performed with 
a period r which is rapid when compared with the char- 
acteristic decoherence time r y~ 1 the probability for more 
than one quantum jump to occur in between two suc- 
cessive projective measurements becomes negligible since 
this is proportional to (jt) 2 , 19 1. The detection and uni- 
tary inversion of single quantum jumps then stabilises the 
quantum information stored in the system. 

The stabilisation scheme presented here allows the de- 
terministic (unitary) restoration of quantum information 
which is to be contrasted with previous proposals where 
the correction contains probabilistic (non-unitary) ele- 
ments PI. 



II. THE TRAPPING MODEL 

Here we briefly introduce the trapping model that we 
use for the analysis of our stabilisation scheme. This 
includes the description of the quantised motion of the 
trapped ion, the relevant internal electronic structure and 
the associated decoherence effects. 





|F = l,mf = 1) 



FIG. 2. (a) Three-level electronic configuration of the 
trapped ion considered here. The two low lying states \a) 
and \b) are coupled through a stimulated Raman interaction 
driven by two (correlated) lasers at frequencies u> ac , Wf, c - The 
upper state |c) is detuned by a large enough amount A that 
it is never appreciably excited, (b) Energy levels of a Be + 
ion which realizes the effective three-level configuration de- 
picted in (a). The states \a) and \b) are associated with the 
\F = 2, m? = 2) and \F = l,rriF = 1) states within the 
fine-structure multiplet of the 2s 2 S\/2 ground state of Be. 
The excited state \c) is associated with the Zeeman-sublevel 
\F = 2, m F = 2) of the 2p 2 P 1/2 excited state [8]. 

Consider the quantized CM motion of a single ion 
which is confined within a trap potential that can be 
closely approximated as being harmonic and which is 
characterized by the three frequencies of oscillation, v Xl 
v y , and v z along the x, y, and z directions fl30| , |3l|| . As 
mentioned above we propose to encode the qubit in the 
two bosonic degrees of freedom which are represented by 
the ion's harmonic motion along the x and y directions. 
We restrict our analysis to those two motional degrees 
of freedom since the ion's motion along the z direction 
plays no role in our scheme. The only decoherence effect 
that we include in our model is motional decoherence. 
To avoid any loss of quantum coherence within the elec- 
tronic degrees of freedom we require a scheme in which 
spontaneous emission is highly suppressed. We therefore 
assume that the relevant electronic degrees of freedom 
form a A system as depicted in Fig.||(a), and that the 
electronic dynamics are confined to the two states |o) and 
| b) which we choose to be two ground state hyperfine lev- 
els. These are coupled by Ml and E2 transitions at best 
so that we can safely neglect any spontaneous emission 
between those states At the same time the \a) \b) 
transition can be driven by two laser beams connecting 
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the ground states \a) and \b) through a common excited 
state |c) in a stimulated Raman scheme [^3[|4|, which we 
briefly review in Appendix [A|. To give a specific exam- 
ple for the realization of an effective three-level scheme as 
indicated in Fig.|^(a) we have drawn in Fig. ||(b) the rele- 
vant energy level structure of Be which, for example, was 
used in a series of experiments in the group of Wineland 
at NIST |,|,@. 

Because the trap potential is harmonic, the position 
and momentum operators can be written in terms of cre- 
ation and annihilation operators for the trap quanta 



2v x m 



Px = i 



hv x 



(&1 -a x ) , 



and similarly for y, and where m is the mass of the ion. 
The free Hamiltonian of the system then takes the form 



j=x,y 



+ s- Sj 



(7) 



where Huj a and Hcob denote the energy of the two ground 
state hyperfine levels and huj c is the energy of the elec- 
tronic excited state. 

We assume that the motional degrees of freedom in the 
x and y directions undergo finite-temperature dissipative 
dynamics which are described by a master equation of 
the form given in Eq. ([l]), and where the associated de- 
cay channels are described by the operators 



1) a T 



C 

x,y 



(8) 



and the quantity n gives the mean excitation number of 
the reservoir Note that for simplicity we have as- 

sumed the decay rate 7 to be the same for the x and y 
directions. In Section [II, we show how a qubit can be en- 



coded in the motional degrees of freedom, such that the 
quantum jumps associated with the decay channel opera- 
tors given in Eq. (||) can be unitarily inverted. In Sections 
IV and ^ where we present our specific proposal for the 
implementation of a stabilisation scheme in a trapped- 
ion system we specialize to the zero-temperature case 
where n = 0, and we are left with the channel opera- 
tors C x _ y = y/jd Xi y. We have chosen this simple form 
of decoherence to exemplify our general strategy and to 
give an explicit implementation for all the steps in our 
stabilisation scheme. 



III. DISSIPATIVE SYSTEM EVOLUTION AND 
STABLE BASIS STATES 

In this section we discuss the dissipative system evo- 
lution and, given the form of the master equation as in 
Eq. (Q) together with the channel operators in Eq. (||), 
deduce the orthogonal basis states \tp±) which allow the 



unitary inversion of the associated decoherence effects. 
We refer to such states as stable. 

Before we go into more detail about the dissipative sys- 
tem evolution, we note that for the derivation of stable 
basis states for a finite-temperature master equation it 
is sufficient to consider the zero-temperature case. The 
reason for this is that the unitary reversibility of a quan- 
tum jump associated with the decay channel C x implies 
the unitary reversibility of a quantum jump through the 
decay channel C' x . This becomes clear from the follow- 
ing argument. If a subspace H.q satisfies the condition 
given in Eq. (||) for the quantum operation A — C x which 
describes a quantum jump through the decay channel 
C x , then the subspace Ho also satisfies the condition in 
Eq. (|J) for the quantum operation A = C' x associated 
with the decay channel C' x . This follows directly from 
the form of the channel operators as in Eq. (||) and the 
commutation relation [a x , a x ] = I, for the bosonic cre- 
ation and annihilation operators. The same argument 
holds for the decay channels C y and C' y . In the following 
we will therefore only consider the zero-temperature case 
where the dissipative system dynamics are described by 
the master equation 



dt • 



■E 

j=x,y 



l -[c]C 3 p + pC}C j -2C j pC}} , (9) 



with the channel operators C x . y — y/ja Xiy , and where 
we have transformed into the interaction picture of Hq 
which is given in Eq. (Q) . For simplicity we have further 
assumed that the system is only subject to its free dy- 
namics. 

We are interested in the decoherence effects during the 
time interval t G [to , to + t] , in between two of the suc- 
cessive projective measurements which serve to detect 
quantum jumps (see Fig.[l]). To examine these decoher- 
ence effects in more detail consider the system evolution 
starting from the initially pure state p(to) = \ip e i) {ip e i \ ® 
where \ip e i) denotes the initial electronic 
state which we include for completeness here although 
it remains unaffected by the dissipative evolution. The 
initial vibrational state \ipmb) = C-\ip-} + c+\ip+) encodes 
a qubit in the two basis states |^_) and which span 
the subspace Tio C TC, of the ions motional Hilbert space. 
During the time interval t £ [to, to+r], we do not monitor 
the dissipative system evolution and the system dynamics 
suffers motional decoherence as described by the master 
equation in Eq. (^). Given this master equation the time 
evolved density operator p (t) can be expressed with the 
help of a Dyson expansion ji4l , [35| 

H*) = E E 

»=° {jija,...,M 

x \t){t - t n )C Jn ...U(t 2 - H)6 h ti{ti - ^ 
p(toP i (ti-t )Cl...W(t-t n )} , (10) 





/•*3 rt 


/ dt n .. 


/ dt 2 \ 


n 


Jt Jt 
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where we sum over all possible sequences of {j n , j±} 
and where ji € {x,y}. In this expression the density op- 
erator p (t) is written as a mixture of pure states which 
result from the non-Hcrmitian time evolution 



U(t) = cxp 



- H H eff t 



(11) 



governed by H e g = —ih(C x C x + C^C y )/2, and inter- 
rupted by quantum jumps through the action of the op- 
erators C x ,y = \fl^x.y We can rewrite Eq. (|Io| ) as a 
sum over sub-ensembles Pn x n^ 

(t) which are character- 
ized by the number of quantum jumps n x , y that these 
sub-ensembles have suffered through the decay channels 
C x ^ y in the time t without specifying the exact times of 
the jumps p5[, 



oo 

= Pn*n y \lpn m ny(t)){l/>n x ny(t)\ 
®\4>el){i>el\- 



(12) 



Here the states \ip nx n [t)) are normalized, so that the 
quantities p Mx n B give the probabilities for the state 
\ipn^n y (t)) in the mixture. To leading order in yt we 
havep„ x „ cx {^t) n:c+nv , pq| . We keep the sub-ensemble 
decomposition as in Eq. (p.2[) and for times jt <C 1 we 



neglect all p„ x „ of order 



and higher to obtain 



P ( T ) = Po x o y P o x o B (r) + Puo y P i x o H (t) 

+Po x i h /5o x i b (t), (13) 

where we have let to = 0, for convenience. 



The mixture in Eq. (13) is of the form given in Eq. 



and serves as our starting point to deduce the stable or- 
thogonal basis states \tp±}- The three terms in the mix- 
ture have the following interpretation. The first term 
corresponds to the system evolving without any quantum 
jumps, whereas the second and the third term describe 
the dissipative effects of a single quantum jump through 
the decay channel C x and C y , respectively, without speci- 
fying the exact time of the jump. From Eq. (|Io| ) the three 
contributions to the mixture can be further evaluated to 
give 



VPiXlV'i.oJr)) 



e< T - 1 



eT r - 1 



C x U(T)\ip vib ) 



CyU(T)\^ mb ) 



(14) 



(15) 



(16) 



where \ipo x ,o y (r)) € Wo x ,o y , \i>i x ,o v (r)) € Hi x ,o y , 
IV'Ox.iy^)) e Kox.iy, and U(t) has been defined in 



Eq. (|ll|). Given that these three subspaces of the mo- 
tional Hilbert space are mutually orthogonal, the mixture 
in Eq. jl^ ) can be projected into one of the three terms 
as we will show in Section IV below. We now consider 



the conditions for the unitary reversibility of the decoher- 
ence effects associated with the three possible outcomes 
of this projective measurement as described by Eqs. jl^ ) 
- ©. 



A. Reversibility of No-Jump Evolution 

First, consider the dissipative system evolution in the 
absence of quantum jumps as described by Eq. (^). Fol- 
lowing Eq. (^) the associated decoherence effects are uni- 
tarily reversible if 



V|V>> e Wo- 



(17) 



and where H.o is spanned by \tp±}- It is not possible to 
find two orthogonal states \ip±), which satisfy this con- 
dition for arbitrary times r, if they are restricted to a 
single degree of freedom of the ion's motion. This be- 
comes clear from the following argument. Let the two 
basis states \<ip±) — J2n c n\ n ): where \n) denote Fock 
states of the ion's motion along, say, the x direction. For 
the decoherence effects in the absence of quantum jumps 
to be reversible for those basis states we will require from 
Eq.(HT 



± e 



which implies |c+| = |c~|, Vn, and the state \ip-} can 
be written as \ip-) — J2 n c n e^"|n). Now, for the basis 
states \ip±) to be orthogonal the phases <f> n cannot all 
be equal and there exists a pair of phases <p ni ^ <f> n2 . 
We can then define a superposition state — — 
e -^"i / \/2, which obviously satisfies € Ho but 
for which Eq. ([Tt]) is not satisfied as can be easily shown. 
We have performed a numerical search and found orthog- 
onal single-mode states which satisfy Eq. ( |l7| ) for specific 
times r. However, these states depend on the particular 
value of t and 7 which limits their applicability and we 
will not present any further analysis of such states here. 

For the orthogonal basis states to be stable and inde- 
pendent of the period r of the projective measurements 
that serve to detect quantum jumps, they have to com- 
prise of at least two bosonic modes. To derive such states 
in two of the ion's motional degrees of freedom we impose 



(18) 



on the basis states \ip±), which span Hq. This choice sat- 
isfies Eq. (17) and ensures that the dissipative evolution 
in the absence of quantum jumps is unitarily reversible. 
In fact, Eq. (03) implies 



G 



(19) 



so that the subspace Ho remains invariant and Hp ,p 
Hp. The general form for states satisfying Eq. Jlg| ) 

\n~, 



• = ^2n=o °n \ n i N — n), where T = jN/2, and where 
n y ) denote the usual number state basis for the ion's 



is 



harmonic motion in the x and y directions. 



are Fock states of the transformed operators a g — (a x + 
e^cLy)/ y/2 and 0,4 = (a y — e^^a x )/ V2, which correspond 
to the transformation of Eq. ( |2l| ) for 9 = tt/4. In the 
Cartesian number state basis the states \tp±) are given 
by 



l^±> = ^{|2,0) 



e- 2 ^|0, 2) ± V2e~^|l 



,!>} 



B. Reversibility of Quantum Jumps 

To continue our analysis of the stable basis states we 
consider the conditions for the unitary reversibility of 
the decoherence effects associated with a single quan- 
tum jump as described by Eqs. and ([if]). Following 
Eq. (Q) these decoherence effects are unitarily reversible 
if 



IVO = A,y 



V|V> e n 



(20) 



where we have used Eq. (|19|). These conditions were sat- 
isfied in the CQED situation of Ref. jlTj by choosing 
= |2,0), and \ipJ) = |0,2), together with the fact 
that the master equation, as given in Eq. (^|) , is form in- 
variant under the general Bogolyubov transformation 



Cx 



V7 



e** sin^ 



(21) 



and the operators C x , v fulfill the conditions of Eq. ( |20| ) 
for 9 = 7r/4. In Ref. fll7[ ] the detected quantum jumps 
were associated with the operators C Xiy , for the values 
9 = tt/4 and = 0. In general, for the same underly- 
ing master equation, the quantum trajectories formalism 
yields different unravelings for different forms of the op- 
erators C x .y The particular values of the parameters 
9 and <fi are determined by the associated measurement 
strategy which monitors the dissipative system evolution. 
For 9 = 7r/4, these measurements yield no information on 
the superposition of the basis states \i/j±), and, so long as 
one can physically realize such stabilizing measurements, 
one can unitarily invert a quantum jump which occurs 
though either the C x or C y decay channel. In the scheme 
of Ref. [0, the measurement which stabilizes the un- 
known quantum superposition against a single quantum 
jump is accomplished through the use of a beam-splitter 
|36f . In our case, however, to perform an analogous mea- 
surement, that of detecting the loss of rotational quanta 
(when 9 = ir/4 and <f> — 7r /2), is difficult. Instead we 
make use of a duality which exists between the form of 
the basis states and the stabilizing measurement. We 
choose, 

|V+> = |0,2) rot , |^_> = |2,0}„ t , 
where the states 

\n d ,n g } ro t = —== a\ a\ |0,0), 
\jn g \n d \ 



For cj) = 7r/2, these are Schwinger two-mode rotational 
states. One can show that for C x _ y — yf^a, x ,y the above 
basis states \ip±) satisfy the conditions given in Eqs. (|17|) 
and (^0|). However, as we have discussed in Section [Tb|, 
to accomplish the detection of quantum jumps we re- 
quire that the decoherence effects described by Eqs. ( |l4| ) 
- (|l6|), must lead to mutually orthogonal subspaces of 
the motional Hilbert space. For the above basis states 
\ip±) this is not the case since they lead to H\ x .o v = 
Ho x ,i = span {|1, 0), |0, 1)} . The orthogonality require- 
ment can however be satisfied by doubling all the excita- 
tion numbers in the basis states, and we define 

|^±)o = ^{|4,0)+ e *|0,4)±V2e*|2,2)} , (22) 

where we have given the most general form with non- 
zero relative phases </>i and 02- These states are sta- 
ble and represent the main result of this section. They 
span a subspace Ho which satisfies the conditions given 
in Eqs. (|l7|) and (p(i|), so that the decoherence effects 
associated with either a single jump or no jump, are 
unitarily reversible. Given an initial vibrational state 
iV'-uife) — C-\i/j-)o + c+|?/>+)o which encodes an arbitrary 
superposition in the stable basis states \tp±)o, we find 
from Eqs. (@) - @ 



c-\ip-) x + c + \i>+) x , 



l^uCr)) = c-\ip-) y + c+\ip+) 



+/y ' 



(23) 
(24) 
(25) 



and the probabilities in the mixture of Eq. (|3|) can be 
evaluated to give p 0:c o y = e" 47T and p 1;r0y = = 
2 e -47T(* e 7T _ 2) decoherence effects lead to mutu- 

ally orthogonal subspaces which are given by Hp x p v = 
span{\ip + )o,\ip-) }, Hi x0y = span {\ip + ) x , \ip-) x } and 
Hp sc \ y = span {\ip+) y , \ip-} y }, where we have defined 



1 



|V±)« = -^{|3,0)±|1,2)}, 
|V±)» = ^{|0,3)±|2,1)}, 



(26) 
(27) 



and we have taken . 



'1.2 



in the definition of the stable 



basis states in Eq. (22). Due to the very special proper- 
ties of these basis states, the decoherence effects cause the 
unknown superposition state \ipvib) — c -|^-)o + c+|VM-)o 
to evolve into a mixture of superpositions with identical 
coefficients c±, but in orthogonal subspaces. If, through 
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the proper measurement, we learn only which Hilbert 
subspace the system resides, then that subspace possesses 
an exact copy of the original qubit state and through a 
unitary manipulation we can restore this information to 
Hq. This does not contradict the no-cloning theorem p7\ , 
since once we have made the measurement, the "copies" 
of the qubit in the other Hilbert subspaces are lost. 

We emphasize that the basis states derived in Eq. (|22 
are stable in the presence of a finite-temperature reser- 
voir even though for their derivation it was sufficient to 
consider the zero-temperature case. It is straightforward 
to show that the subspace Ho = span {\ip+)p, \ip—)o}, also 
satisfies the conditions given in Eqs. (O) and (EOT) for the 
finite-temperature case where the associated decay chan- 
nel operators are given by Eq. (|^) . In this context we note 
that recent work on the more formal aspects of states of 
bosonic systems that can be stabilized in the presence of a 
zero-temperature reservoir has also uncovered the states 
given in Eq. (|22| ) and many more complicated multi-mode 
states Q . In addition to the results reported in Ref . |3|] 
the analysis presented here demonstrates the potential of 
those states for the more general, finite-temperature case. 
In the remainder of this paper we concentrate on the case 
of a zero-temperature reservoir as we have already noted 
above. 



IV. DETECTION OF DECOHERENCE 
PROCESSES 

In the following we discuss how to detect the decoher- 
ence processes associated with the master equation given 
in Eq. (|9|). More specifically, based on our specific choice 
of encoding information in the motional basis states de- 
rived in Eq. (j22]), we present a method that determines 
whether a quantum jump associated with the decay chan- 
nel C x or C y has occurred. The corresponding processes 
constitute step (i) in our stabilisation scheme as depicted 
inFig.0. 

The detection of a quantum jump is accomplished in 
two steps. The first step serves to generate an "error 
syndrome" which only carries information about the ef- 
fects of decoherence. To be more specific, given that the 
initial e lectro nic state of the ion \ip e i) — \ a) , we show in 
Section IV A how to entangle the mixed state that re- 
sults from the decoherence processes, as in Eq. (|l3|), with 
the electronic states \a) and \b). In essence, we generate 
a unitary operator U en t, which, in its generic form, is 
described by 



Uentllpvib) <S> | O) = \lpvib) <S> \b) , y\lpvtb) € Hent , 
Uent\lpvib) <& | a) = IV'mft) ® I a) , V\lpvib) € U^ nt , 

(28) 

where H en t is a specific subspace of the ion's motional 
Hilbert space and Ttjr nt is the orthocomplement of Hent- 
It is important that this operation leaves the motional 



state \ip v ib) unaffected. Also it is necessary that the op- 
erator U en t entangles the motional subspaces Ttjr nt and 
TL eri t with the electronic states \a) and \b) without sup- 
plying any further information on the vibrational state 
\ipmb). This can be accomplished by making use of our a 
priori knowledge of the motional basis states and the de- 
coherence processes discussed in Section III. The above 



entanglement process produces a binary signature in form 
of the ion's electronic state which serves as our error syn- 
drome. The second step of our quantum jump detection 
scheme serves to "read out" this error syndrome. This is 
accomplished through a projective measurement of the 
ion's electronic state which we discuss in Section |IV B| . 
Following the entanglement operation given in Eq. (|28]), 
the detection of the ion in the electronic state \a) (\b)) 
amounts to projecting the motional state of the ion onto 

[Hent)- 



the subspace Hjr nt 



Before we describe in detail how the error syndrome 
can be generated and read out we briefly elaborate fur- 
ther on the sequence of operations which constitute step 
(i) in Fig.^. First, the detection of a quantum jump 
associated with the decay channel C x is accomplished 
through constructing the entanglement operation such 
that Hent — "^l^Oy in Eq. (pgj). A subsequent measure- 
ment of the ion in the electronic state \ip e i) — \b) is then 
equivalent to the detection of a quantum jump as de- 
scribed by Eq. (plf), and the motion of the ion is pro- 
jected into the pure state given in Eq. (24). If, on the 
other hand, the projective measurement yields the re- 
sult \ipei) = |o)j for the electronic state of the ion, no 
quantum jump associated with the decay channel C x has 
occurred and the ion is left in a mixture of the states 
given in Eqs. (|23|) and (|25|). Irrespective of the out- 
come of this first projective measurement we then de- 
termine whether or not a quantum jump associated with 
the decay channel C y has occured. The corresponding- 
procedure is identical to the one which serves to de- 
tect a quantum jump through the decay channel C x but 
where now Hent — HQ x i y for generating the error syn- 
drome. If through the first measurement we have de- 
tected a quantum jump associated with the decay chan- 
nel C x the outcome of the second measurement will al- 
ways be \ip e i) = | a) and we do not obtain further infor- 
mation. However, for the case where no quantum jump 
through the decay channel C x has occured, the outcome 
\ipel) — \b) for the second projective measurement corre- 
sponds to the detection of a quantum jump as described 
by Eq. (^|) and projects the ion into the pure motional 
state given in Eq. ( j25[) . On the other hand the measure- 
ment of the ion in the electronic state \ip e i) — |a) for a 
second time implies that no quantum jump has occurred 
through either of the decay channels and the ion is pro- 
jected into the pure motional state given in Eq. (|2^) . The 
second projective measurement thereby completes the se- 
quence of operations that constitute step (i) in Fig.[j] and 
leaves the ion in a definite, pure motional state. From 
this the original motional state can be restored in step 



8 



(ii) of our stabilisation scheme as indicated in Fig. [j] and 
as we describe in Section 0. In the following we now de- 
tail the processes which accomplish the generation and 
reading out of the error syndrome. 



A. Entanglement Stage 

The physical process that we use to generate the essen- 
tial part of the entanglement operation given in Eq. fl2S| ) 
is a resonant stimulated Raman transition between the 
two ground states \a) and \b) in a specific laser geometry. 
This generates Rabi oscillations between the states \a) 
and | b) where the Rabi frequency depends on the ion's 
motional state. An essential requirement for the genera- 
tion of the desired entanglement is the accurate control of 
the phases acquired through these Rabi oscillations. We 
therefore consider time dependent laser pulses character- 
ized through their dimensionlcss pulse shape f(t) to gen- 
erate the Raman transitions, and the phase of the Rabi 
oscillations is controlled through the generalized pulse 
area 



.4 



f(t'fdt'. 



(29) 



For the generation of the error syndrome associated with 
the decay channels C x and C y , we will be interested in the 
case where the Raman-induced Rabi oscillations are only 
sensitive to the ion's motion along one of the principal 
axes. As shown in Appendix A|, this can be realized by 
resonantly exciting the two-photon transition \a) \b) 
with two plane waves having wavevectors fc a and k b , and 



where the wavevector difference 5k 



k b is aligned 



with that axis. The resulting Hamiltonian is given by 
-%(i)exp [-tfj/2] \a)(b\ 



H 



OO 

n=0 



,!nl 



H.c. 



(30) 



where j = x for the case of aligning 5_k with the x axis, 
and j — y for the case of aligning 5_k with the y axis. The 
symbols rj x and rj y denote the Lamb-Dicke parameters 



Tlx = Ax 



(31) 



for the case j = x, y, as explained above, and Axo = 
(H/2v x m) 1 / 2 and Ayo = {h/Zvym) 1 / 2 give the widths of 
the motional ground state along the x and the y direction, 
respectively. The Raman coupling constant g(t) oc f{t) 2 
and its exact form is given in Eq. (A4). To gener- 
ate the required entanglement operations we assume the 
Lamb-Dicke limit [P7lf28|, which can be defined through 



(ip vlb \5k 2 x 2 \i! 




and (ip V ib\5kyy 2 \ip vi b) < 1, and 



which implies rj x ,y 1, for the Lamb-Dicke parameters. 
This allows us to retain only the lowest order terms in 
r\ x . y in the Hamiltonian describing the Raman-induced 
dynamics. 



We first address the question of detecting whether a 
quantum jump associated with the operator C x = y/^a x 
has occurred. The corresponding error syndrome is gen- 
erated through an entanglement operation as in Eq. (|2q) 
with H en t = Hi x o y , and H^ nt = Wo x o y © W 0x i H , which 
can be accomplished by exploiting the fact that for all 
states \ipvib) € 7ii x o y the motional excitation numbers 
n x along the x direction are odd, whereas for states 
\ijj v ib) € Ho x Oj, © 7^o x i B the motional excitation numbers 
n x are even, as seen from Eqs. (p2|) , ( p6| ) and (f27|). With 
this in mind, we construct a laser exciting field where the 
induced stimulated Raman transition is only sensitive to 
the ion's motion along the x direction. More specifically, 
we realize the essential part of the above entanglement 
operation through the unitary time evolution generated 
by the Hamiltonian 



H x = hg(t) a x a x <g> \a)(b\ + H.c. 



(32) 



for a specific choice of the generalized pulse area given in 
Eq. (|29|) . This form of Hamiltonian has previously been 
discussed as a degenerate Raman-coupled model in the 
context of CQED (39). More recently, Gerry has shown 
how such an interaction may be realized for a trapped 
ion by driving a dipole transition in a specific laser ar- 
rangement pofl. 





JJ R (2) R (2) 
5kW 
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R 1 




(a) (b) 
FIG. 3. Laser arrangement which generates the error syn- 
drome for a quantum jump associated with the decay chan- 
nel C x - (a) Two pairs of lasers drive a stimulated Raman 
transition between the states \a) and |6), which is sensi- 
tive to the ion's motion along the x direction. The first 
pair, is arranged so that the wavevector difference 

5k^ — kS^ — k^\ is aligned with the x axis. The second pair 
R i 2 b is co-propagating, so that 5k (2) = k^ } - fcf > = 0. (b) The 
stimulated absorption and emission processes induced by the 
laser excitation can be treated separately for the two pairs 
of Raman lasers, if the detunings A^ 1 - 1 and A*- 2 ' are chosen 
sufficiently different. 

For the dipole-forbidden \a) <^> \b) transition which is 
relevant to us here, the form of coupling given in Eq. (p^) 
can be generated through a resonant stimulated Raman 
transition between the electronic states \a) and \b) which 
is induced by two pairs of laser beams as shown in Fig.|^. 
The first pair of Raman lasers is arranged so that the 
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wavevector difference is aligned with 

j) this generates the Hamil- 



the x axis. Following Eq. 
tonian 



HW = -hg^(t)exp [-jf/2] \a)(b\ 
(g) { 1 - T) 2 a\ a x } + H.c. . 



(33) 



where we have assumed the Lamb-Dicke limit [ p7| , and 
kept only the leading order terms in the Lamb-Dicke pa- 
rameter 77 = rj^K The second pair of Raman lasers serves 
to cancel the zeroth-order contribution in this Hamilto- 
nian which is insensitive to the motional state of the ion. 
The laser geometry is arranged such that the wavevector 



difference Sk^ = k^ 2 



..(2) 



0, vanishes exactly |tl| 



From Eq. (|30|) this gives rise to the Hamiltonian 



ffW = -hg^{t)\a)(b\ +H.c. 



(34) 



(2) 

since rjx — 



(2) 

Vv 



= 0. The combination of the two 
pairs of Raman lasers is described by the Hamiltonian 
H = + H^ 2 \ if the stimulated absorption and emis- 
sion processes induced by the laser excitations can be 
treated separately for the two pairs of Raman beams. 
This is the case if the detunings A^ and A^ 2 ) of the 
lasers from the excited state |c) (Eq. (A2)) are chosen 
sufficiently different for the two pairs of Raman beams 
as depicted in Fig. ||(b). In this situation, and if we as- 
sume the laser phases and amplitudes to be arranged 
such that g^ 2 \t) = — exp [— rj 2 /2]g^ (t), the laser excit- 
ing field described above generates the Hamiltonian given 
in Eq. (|32|) and the coupling constant 



g(t) = V 2 e X p{-r, 2 /2}g^(t)=gf(t) 2 



(35) 



For the Hamiltonian in Eq. (|3^) one can calculate the 
unitary operator which describes the resulting time evo- 
lution to be 



U(A)*=c 8\xA](\a)(a\+\b)(b\) 

+ sm[xA}(\a)(b\-\b)(a\) 



(36) 



where \ = |<7| o-J. a.^, and we have further assumed 
g = i\g\ for the phase of the coupling constant. The 
generalized pulse area A is given by Eq. ( p9| ) . If we now 
choose A = n/2\g\, we obtain the mapping 



U(A)*\1>±) X <8> \a) = KM*®|6>, 



(37) 



for the basis states \if>±} x which are given in Eq. ( |26| ) 
and which span the motional subspace 7ii x o H associ- 
ated with the decay channel C x . On the other hand, 
for the basis states \ip±)o and \ip±) y which are given in 
Eqs. ( |2^ ) and (pTj), and which span the motional subspace 
7io x Oj, ©Hq x i , the time evolution generates the mapping 



U{A) x \ijj±}o iy ® \a) = \ife)o, v ® |o) 



(38) 



so that the laser arrangement described above generates 
the essential feature of the entanglement operation given 



in Eq. ([28|) . It transfers the binary information as to 
whether or not a quantum jump associated with the de- 
cay channel C x has occurred, into the electronic degrees 
of freedom. However, the Raman-induced operation does 
not leave the motional state of the ion unaffected. As 
seen from Eqs. (^) and ([58]), it interchanges the basis 
states that encode the quantum information in the sub- 
spaces Ho^Oy j H-i^Oy and 7~to x i v ■ However, after determin- 
ing whether or not a quantum jump has occurred through 
the C x channel, we must repeat the interrogation to de- 
termine whether or not a jump has occurred through the 
C y channel. This second interrogation again interchanges 
the motional basis states and thus after both steps are 
completed, the motional state is returned to its original 
configuration. 



R;, 



R (l) 



/ 



V R (2) R (2) 
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(a) (b) 
FIG. 4. Laser arrangement which generates the error syn- 
drome for a quantum jump associated with the decay channel 
C y . (a) The laser geometry is almost identical to the one de- 
picted in Fig. ^ but where the first pair R^j is now arranged 

so that the wavevector difference Sk} 1 ^ = fe^ 1 ' — k^ , is aligned 
with the y axis, (b) The laser frequencies remain the same as 
in generating the error syndrome for the C x decay channel. 



The detection of the quantum jump associated with the 
decay channel C y is achieved through an almost identical 
procedure. The corresponding error syndrome is gener- 
ated through the entanglement operation 



U{A) y \iP±) , x ® \a) = \ip T ) ,x <8> I a) 
U(A)y\^ ± ) v ® \a) = hM»®|&>. 



(39) 



This mapping is generated in an almost identical manner 
as the one which serves to generate the error syndrome for 
the decay channel C x . The unitary operator U(A) y can 
be realized through resonantly driving the \a) \b) tran- 
sition with two pairs of Raman lasers as described above 
but in a slightly different laser geometry. This is illus- 
trated in Fig.[|. The first pair of lasers is now arranged 
so that the wavevector difference SkS 1 ^ 1 = k^ — k^\ is 
aligned with the y axis, and the resulting Raman tran- 
sition is only sensitive to the motion of the ion along 
the y axis. The geometry of the second pair of Raman 
lasers and the arrangement of the laser phases and am- 
plitudes remain the same, so that with the replacement 
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x — > y the analysis is identical to the discussion that led 
us from Eq. ( |33| ) to Eq. ( ^6|) . The resulting time evolution 
is described by 



U(A)« = co S \xA](\a)(a\ + \b)(b\) 

+ sm[xA](\a)(b\-\b)(a\) 



(40) 



where \ — Iff I &y an d the coupling constant g is de- 
fined in Eq. (pq). With the choice A = 7r/2|g|, the time 
evolution generates the mapping given in Eq. ([59]) and 
the error syndrome for the detection of quantum jumps 
associated with the decay channel C v can be read out 
through a measurement of the ion's electronic states \a) 
and 1 6) as we will describe in Section IV B below. 

To give an estimate of the time scales involved in the 
entanglement stage discussed here, we consider the spe- 
cific example of Be + ions and associate the electronic 
states \a), \b) and |c) as indicated in Fig.||(b). This is 
the level scheme that was employed by Monroe et al. , 
for the first demonstration of quantum logic with trapped 
ions and in further experiments in the group of Wineland 
at NIST Ji|,[l0| . To remain specific we model the explicit 
time dependence of the laser pulses through the dimen- 
sionless pulse shape 



f(t) = sin 2 



irt 
Tl 



0<t<T L , 



(41) 



where Tl denotes the pulse duration. We can then 
evaluate the generalized pulse area in Eq. ( p9] ) to ob- 
tain A = 3Tl/8. The generation of the error syndrome 
through the entanglement operations in Eqs. - ( pij ) 
requires A = iv/2\g\, for the pulse area which implies 
Tl = Aiv/3\g\, for the pulse duration. To estimate this 
time we calculate the value of \g\ from Eq. (|35|) . We 
assume \g^\/2n = 500kHz for the peak value of the 
coupling constant generated by the first pair of Raman 
lasers as in Eq. ( |33| ) and n = 0.2, for the correspond- 
ing Lamb-Dicke parameter. These are parameters taken 
from the experiment reported in | fL0| . From Eq. (|35| ) we 
find |<?|/2tt = 20kHz, and obtain 



T L ~ 33^s , 



(42) 



for the duration of the pulse which generates the error 
syndrome. This can be shortened through an increase in 
the value of \g\. However, to remain within the validity 
of our analysis we have to consider the limit set by the 
vibrational rotating wave approximation that was made 
in Appendix ^| in deriving the Hamiltonian in Eq. (30). 
This limit is given by the off-resonant excitation of the 
first vibrational sideband through the first pair of Ra- 
man lasers and requires rj\g^ \/u x>v <C 1, ]42j] . For the 
quoted experimental parameters and u x _ y /2iT — 10MHz, 
jUjl , we find "q\g^ \/v x ,y — 10 -2 , so that the duration of 
the entanglement operation may be shortened by one or- 
der of magnitude. To go further would however require 
an increase in the trap frequency to remain within the 
low-excitation regime. 



It is important to contrast the entanglement oper- 
ation described here with the transfer of a quantum 
state between atomic ground-state Zeeman sublevels and 
bosonic degrees of freedom as described by Parkins et al. 
through adiabatic processes in a CQED situation p3[ . 
There, a quantum state is transferred, whereas the entan- 
glement operation described here leaves the vibrational 
state \ipvib) unchanged and transfers information as to 
whether this states lies in a specific subspace of the mo- 
tional Hilbert space. 



B. Projective Measurement 

We now turn to the second step of our quantum jump 
detection scheme which reads out the error syndrome 
generated through the entanglement operation discussed 
above. As we have already mentioned, this second step 
consists of a projective measurement of the ion's elec- 
tronic state. This measurement must not significantly 
disturb the motional state of the ion which we aim to 
stabilise here. Existing measurement procedures such as 
atomic shelving and resonance fluorescence [ p3|]24| ] entail 
the emission of a large number of spontaneous photons 
in the typical experimental situation where detector effi- 
ciency is less than unity and the solid angle covered is not 
the full 47r. In these measurement strategies the center of 
mass of the ion receives numerous spontaneous emission 
recoil kicks in random directions which results in an over- 
all heating of the ion's motion. To avoid this we employ 
a technique which has been discussed in the literature for 
the generation of single photon wavepackets |2^,Q . We 
adopt here the scheme of Gheri et al. ^6|, where a A- 
type three-level atom (in our case the ion whose internal 
level structure is shown in Fig.||) is surrounded by a low- 
Q optical cavity and excited by a time-dependent laser 
pulse. The laser and cavity electric fields couple the two 
ground states \a) and \b) through a resonant two-photon 
process and induce a stimulated Raman transition be- 
tween those states. The excited state |c) is detuned by 
a large enough amount A that it is never appreciably 
excited and spontaneous emission events can be faith- 
fully neglected. The situation is illustrated in Fig.|5|. We 
show below how, with suitable approximations, the ion's 
electronic and motional degrees of freedom decouple so 
that the ion-cavity system reduces to the model consid- 
ered by Gheri et al.. We will not go into the detailed 
solutions to the dissipative dynamics since these are dis- 
cussed in Ref. (^6|. Gheri et al. have shown that if the 
initial electronic state \ip e i) = \b), then the laser pulse 
transfers the electronic population into state \a) and si- 
multaneously excites the cavity mode which rapidly de- 
cays into the external field in the form of a single photon 
wavepacket. More specifically, in the bad cavity regime 
where k 3> |<?(i)|, they derive an explicit solution to the 
dynamics of the ion-cavity system [^6j . Using this solu- 
tion the probability P(t) for a single photon to be trans- 



it 



mitted from the cavity during the time interval [0, t] can 
be calculated as shown in Ref. Hi and we obtain 



P(t) = 1 - exp 



(43) 



where k is the decay rate of the cavity field amplitude, 
and g(t) denotes the Raman coupling constant for the 
transition between the two ground states \a) and \b). 
Note that g(t) oc E L (t), where E L (t) = f(t)E L denotes 
the electric field amplitude of the external laser pulse and 
determines the time dependence of g(t) through the pulse 
envelope f(t). 



K 





FIG. 5. To perform a measurement of the ion's electronic 
states | a) and |6) without significantly disturbing the mo- 
tional degrees of freedom the ion is surrounded by a low-Q 
optical cavity and excited by a laser pulse with electric field 
envelope E(t) as indicated in (a). The laser and cavity elec- 
tric fields resonantly couple the two states \a) and |6) and 
induce a stimulated Raman transition between those states 
as shown in (b). In the bad cavity regime this transition 
is overdamped. The excited state \c) is detuned by a large 
enough amount A that it is never appreciably excited. If the 
electronic state of the ion \ifr e l) = \b), then the laser pulse 
transfers the population into state \a) and simultaneously ex- 
cites the cavity mode which rapidly decays into the external 
field in the form of a single photon wavepacket. On the other 
hand if the electronic state of the ion |V>ei) = \o), no photon is 
generated. In this situation the presence (absence) of a trans- 
mitted photon which is observed through placing a detector 
in the cavity output channel as indicated in (a), constitutes 
a measurement of the ion in its electronic state |6) (\a)). 

For a sufficiently large pulse area, the exponential term 
in Eq. (|i^ ) becomes negligible and P(t) 1. This implies 
that if the initial state of the ion \tp e i) = \b), then the 
laser pulse triggers the transmission of a single-photon 
wavepacket from the cavity with almost unit probability. 
In this situation the detection of a transmitted photon 
through placing a detector in the cavity output channel 
as indicated in Fig.||(a) constitutes, to a very good ap- 
proximation, a measurement of the ion in its electronic 
state | b). On the other hand, if the initial electronic state 
\ipel) = \a), no photon is generated since the frequencies 
of the cavity and laser electric fields are assumed to be 



significantly different from each other, so that they only 
generate transitions associated with their own channels 
as shown in Fig.||(b). Therefore, the absence of a de- 
tection event represents a measurement of the ion in its 
electronic state \a). Note that irrespective of the mea- 
surement outcome the ion is left in its electronic state 
\ipel) = \ a ) after the measurement. 

In the following we show how the ion-cavity system 
depicted in Fig.|3| reduces, in the low-excitation regime 
J42[ and the Lamb-Dickc limit p7| , ^8[ , to the model con- 
sidered by Gheri et al. [^6| . The Hamiltonian which de- 
scribes the free evolution of the electronic, motional and 
cavity degrees of freedom is given by 



E 

j=a,b,c 



hv i 

o=x,y 



1 



hue A? A, 



(44) 



where the operators A (A T ) are the bosonic annihilation 
(creation) operators for the single-mode cavity field ex- 
citations and loc denotes the cavity-mode frequency. As 
depicted in Fig. ||(b) we are considering a resonant two- 
photon transition between the states |o) and \b) and as- 
sume the triggering laser pulse and the cavity mode to 
be detuned by the same amount A from the excited state 
\c). The corresponding electric fields are given by 



E L (x,t) = E L (t) e -*t*i*-"i*] + H.c. 
E c (y) = E c cos [key] + H.c. , 



(45) 



where ki, denotes the wavevector of the laser and kc de- 
scribes the mode function of the cavity. As seen from the 
expression for the cavity electric field, we assume the ion 
to be positioned at an antinode of the cavity field mode 
to minimize the coupling between the ion's electronic and 
motional degrees of freedom which arises from the elec- 
tric field gradients Q. In the rotating frame defined by 
the Hamiltonian H = Hq — HA\c){c\, where Ho is given 
in Eq. (Q) , the coherent dynamics of the system are de- 
scribed by the Hamiltonian 

Hi(t) = hA\c)(c\ - hg L (t) \b)(c\ e-^Wl + H.c. 

-hg c \a) (c| i f cos [k c y(t)] + H.c. , (46) 

where the coupling constants <7z,(i) = {b\p\c) E ^(t) / h, 
and gc — {a\p\c)Ec/h, in terms of the laser pulse am- 
plitude EL(t) and the electric field per photon Ec inside 
the cavity. The symbol p denotes the dipole operator and 
we have made the optical rotating wave approximation 
in the above. Note that in the rotating frame the po- 
sition operators are explicitcly time dependent and are 
given as x(t) — Axo(a :E e _Wxt + a\, e Wxt ), and y(t) = 
Ay (a y e~ lUyt + a* e Wyt ), where Axq — (h/2v x m) 1 / 2 and 
Ay = (h/2v y m) 1 / 2 , give the widths of the motional 
ground state along the x and y directions, respectively. 
In the limit of large detuning, where |A| ^> |<7l(£)|, |<7c|, 
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we can adiabatically eliminate the excited state |c) 
obtain 



and 



Hj{t) 



-M a \a){a\&A cos 2 [k c y(t)} ~ W) \b)(b\ 
-Hg(t) \a)(b\ it cos [k c y{t)] e ^ £ ^ + H.c. 



(47) 



as shown in Ref. [|34j. The first two terms in this ex- 
pression correspond to optical Stark shifts arising from 
the adiabatic elimination of the excited state and we 
have introduced the abbreviations i5 a = |gc| 2 /A and 
Sb(t) = \gL(t)\ 2 /A. The last two terms describe the stim- 
ulated Raman transition between the two ground states 
\a) and \b), and g(t) — gcg* L {t)//S. denotes the corre- 



sponding Raman coupling constant. To demonstrate how 
the dynamics described by the Hamiltonian in Eq. (|47] ) 
reduce to that of the model of Gheri et al. pq] , we 
first perform the vibrational rotating wave approxima- 
tion [B8[. To this end we expand the position operators 
in Eq7([47|) in terms of the vibrational creation and anni- 
hilation operators a x and a x _ y in normal-ordered form 
Pq| . We then assume the low excitation regime, where 
S a ,Sb(t) <§; v x ,y, and the Stark shifts produced by the 
laser and cavity electric fields are small when compared 
with the vibrational energy level spacing |4^ |. This im- 
plies that the resonance conditions are not significantly 
modified. Keeping only resonant terms in the Hamilto- 
nian in Eq. (W7 ) we obtain 



\2n 



„,„! & l SS\-W)\bM 



Hj(t) = -M a \a)(a\ iU J l + e -*n>a f] 

* — ' n\n\ 

K n=0 

°° f; m \2n /'„•„ \1m 

-hg(t)\a)(b\Ai ® e -^W L ]/2 £ («*?) M 



. . a t ®v + H.c. 



n,m—0 



(48) 



where the only time dependence stems from the enve- 
lope of the external laser pulse. In writing Eq. ([l8]) we 
have further assumed that the corresponding pulse turn- 
on and off is slow when compared with the trap frequen- 
cies v x and v y . This implies that the duration of the 
laser pulse needs to be long when compared with v~y, 
to avoid the excitation of vibrational sidebands through 
frequency components within the bandwidth of the pulse. 
In Eq. (|4§|) we have introduced the Lamb-Dicke param- 
eters r]L = k^Axo and r\c = kcAyo for the transi- 
tions that are coupled through the laser and cavity elec- 
tric fields, respectively. In the Lamb-Dicke limit p7| , p8[ , 
Vl,c ^ lj and we keep only the lowest order terms in the 
Lamb-Dicke parameters to obtain 

fl>(t) = -hS a \a) (a\ A^A - H5 b (t) \b) (b\ 

-hg{t) | a) (b\ i f + H.c, (49) 

so that the motional degrees of freedom decouple and the 
system depicted in Fig. |5| reduces to the scheme of Gheri 
et al. where only the electronic and cavity degrees of 
freedom participate in the dynamics. 

In the bad-cavity regime the cavity decay rate k sets 
the fastest time scale in the system and dominates the 
dynamics of the \a) <=> \b) transition inhibiting any Rabi 
oscillations between the states \a) and \b). In this limit, 
given that the initial state of the ion-cavity system is 
\tp) = |6) |0)c, and the coherent system dynamics are 
described by the Hamiltonian in Eq. (p9|), Gheri et al. 
derive the solution 



\ip(t)) = cxp 



x (|6>|0> c + i2j£|a)|l>c) , 

for the ion-cavity state before the cavity excitation has 
been lost to the external field and where \n)c denote the 
usual number states for the electric field inside the cav- 
ity. With this result the probability of detecting a single 
photon during the time interval [0, t] can be calculated 
as P(t) = 2k J%p{t')\A^ A\tl){t'))dt' , if an ideal photode- 
tector is used p4fl . From this we obtain the expression in 
Eq. (|i^), where the time dependence is essentially con- 
trolled through the triggering laser pulse. 

We have already mentioned that for a near perfect 
projective measurement of the ion's electronic degrees of 
freedom we require 



Arg : 



dt > 1 



(50) 



for the argument of the exponential term in Eq. (|43|). 
This can be achieved through a sufficiently large pulse 
area or pulse duration which we have denoted by Tl, 
here. To further elaborate on this condition we consider 
again the specific example of Bc + ions as indicated in 
Fig. ||(b) . For the value of the cavity decay rate and the 
ion-cavity coupling strength we assume k/2tt = 750kHz 
and \gc\/2T: = 5MHz, which are experimentally realistic 
parameters p6| . To remain specific we model the ex- 
plicit time dependence of the laser-ion coupling through 
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the pulse shape given in Eq. ( |4l| ) so that ffi(t) = f(t)gL, 
and assume \gi\l2-n = 5MHz, for the peak value of the 
corresponding coupling constant. In view of the exper- 
iment reported in Ref. fl(i| this is a moderate assump- 
tion. For the value of the detuning of the laser and 
cavity electric field from the excited state |c), we as- 
sume \A\/2tt = 250MHz > |si,|,|sc|, so that the con- 
dition for the adiabatic elimination of the excited state 
is satisfied. The parameters for \gL,\, \gc\ an d |A| lead to 
\g\/2ir = 100kHz, for the peak value of the Raman cou- 
pling constant and the assumption of the bad-cavity limit 
where \g\ -C k is justified. Furthermore, this value for 
\g\ and the value of the Stark shifts induced by the laser 
pulse and the cavity electric field, |<5 a ,&|/27r = 100kHz, are 
well below the typical trap frequencies v Xty /2ir sa 10MHz, 
in experiments with Be + [To| ], so that the assumption of 
the low-excitation regime is valid. With these parame- 
ters and the expression for the laser pulse shape given in 
Eq. (|il]) the argument of the exponential term in Eq. ( 4^) 
can be calculated explicitely and the condition in Eq. ( 50) 
imposes the limit 



quantum jump associated with the decay channel C x>y 
results in the state 



T L > 



4k 



16/is , 



(51) 



on the duration of the external laser pulse. This condi- 
tion also satisfies the requirement Tl > ~ 16ns, so 
that the bandwith of the pulse is small when compared 
with the separation of the vibrational sidebands. For 
T L = lOOus, the probability in Eq. @, P(T L ) = 99.8%, 
and the detection of the photon transmitted from the 
cavity represents a near perfect projective measurement 
of the ion in its electronic state \tp e i) = \b). Thus, through 
the introduction of a low-Q optical cavity we are able, in 
a particular limit, to read out the error syndrome which 
is stored in the electronic ground states without distur- 
bance of the motional state. The duration of the projec- 
tive measurement can be shortened for larger values of 
the Raman coupling constant \g\, which may be achieved 
through an increase of the electric field of the triggering 
laser pulse, if at the same time the trap frequencies can 
be increased to remain in the low-excitation regime. 



V. THE UNITARY INVERSION 



In Sections [II and IV we have described how to en- 



code the qubit in stable motional basis states and how 
to detect the quantum jumps associated with the master 
equation given in Eq. (0). We now turn to the unitary 
inversion of the detected quantum jumps, the processes 
which constitute step (ii) of our stabilisation scheme as 
depicted in Fig.|l|. In the case where no quantum jump 
has been detected through the interrogation process de- 
scribed in Section |y| the ion is left in the motional state 
given in Eq. which is identical to the initial motional 
state in Eq. (p[), so that no further manipulation is re- 
quired to restore the quantum information to the origi- 
nal basis states. On the other hand, the detection of a 



\lpvib)x,y = C-\tp-) Xj y + C+\lp+) a 



(52) 



corresponding to the channel operators C x and C y , and 
the states \ip±) x ,y have been given in Eqs. ( p6| ) and (p7|). 
The unitary restoration of the quantum information then 
requires the transformation 



Ure V sM±)x.,y 



VP 



±0- 



(53) 



The generation of this transformation is not simple. The 
best solution would be to directly generate the unitary 
operator in Eq. (^) by "dialing up" a specific finite uni- 
tary operator which acts within the two-dimensional mo- 
tional Hilbert space. Although there has been some dis- 
cussion concer ning the generation of any given finite uni- 
tary operator [}48fl , this has yet to be implemented for 
the states of motion of a trapped ion. We were able to 
find a discrete, finite two-mode unitary transformation 
which will effect the restoration in Eq. (53) E9|, but we 



know of no systematic physical mechanism which could 
implement this transformation. In the absence of such a 
mechanism we are forced to look for a particular sequence 
of processes that will effect the unitary inversion of the 
quantum jumps associated with the decay channels C x 



and C y as in Eq. (|53j) . We will attempt to avoid the use 
of resonant interactions as they are prey to severe timing 
constraints. Instead, we will effect the restoration using 
two adiabatic transfer processes and a single, intermedi- 
ate stimulated process. 

We will only describe the procedure which accom- 
plishes the transformation in Eq. ( |53| ) for the case where 
a quantum jump via the C x = ^a x channel has been 

detected. The case for the C y = y^ja y channel is almost 
identical. After the detection of a quantum jump asso- 
ciated with the decay channel C x , the ion is left in the 
product state \ip) x = \ip v ib)x ® |a), where the motional 
state is 

\ip vib ) x = ±={c+ + c_)|3,0) + (c+ - c_)|l,2) , (54) 

and where the coefficients c± are unknown, as seen from 
Eqs. ( ]24] ) and (^6|). The first step in the restoration of 
the state \ip V ib) x to the original state \tp V ib) will be to co- 
herently add one quantum to the vibrational excitation 
in the x direction. To accomplish this we make use of 
adiabatic passage techniques [f50[ , which have previously 
been studied in the context of coherent population trans- 
fer filf . We coherently transfer population between the 
electronic states \a) and \b) through a pair of overlapping, 
time-delayed laser pulses which are known as the pump 
and the Stokes pulse ]5l| . More specifically, we assume 
that the laser which generates the pump pulse is resonant 
with the | a) <£4> |c) transition, and the laser that gener- 
ates the Stokes pulse is tuned to the first red sideband of 
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the 1 6) \c) transition, and aligned with the x axis. As 
shown in Appendix [b|, this excitation of the A system, 
depicted in Fig.||, is described by the Hamiltonian 

Hadd = -hg a (t)\a)(c\ - hg b {t) a) x <g> \b)(c\ + H.c. , (55) 

where the time dependence of the coupling constants 
g a (t) and gb(t) stems from the dimensionless shapes f a (t) 
and fb(t ) of the pump and the Stokes pulse as seen from 
Eq. The above Hamiltonian possesses the particu- 

lar instantaneous eigenstate 



dark 

_ a{t)\n x ,n y ) g \a) - f3(t)\n x + l,n y ) g \b) 



(56) 



whose corresponding eigenvalue vanishes. This eigen- 
state does not contain any contribution from the ex- 
cited state |c), and is generally known as a "dark state" 
|52t . The explicit time dependence of this state is given 
through the quantities a(t) — g^(t)[(n x + l)!/n x !] 1 / 2 and 
(3{t) — g*(t). The addition of a single quantum in the 
vibrational excitation along the x direction is accom- 
plished through adiabatically following the dark state 
in a counter-intuitive pulse sequence, where the Stokes 
pulse overlaps but precedes the pump pulse JHlj ]. In this 
situation the dark state \ipn a n y (*)) dark = \n x ,n y ) (g) \a), 
before the pump pulse is turned on, and \ip nas n v (fy) dark = 
\n x + 1, n y ) ® \b), after the Stokes pulse is turned off, in- 
dependent of the vibrational excitation numbers n x and 
n y . If the pulse sequence is performed adiabatically the 
system follows the dark state given above and, after the 
adiabatic passage, the ion is left in the state 



1 



V2 



{( C+ +c^)|4,0) 



-(c+-c_)|2,2)}® |6) 



(57) 



This is accomplished without severe constraints on the 
pulse durations and amplitudes, as long as the process 
is performed adiabatically. The adiabaticity condition is 
satisfied if 



\9a,b(t)\dt^> 1 



(58) 



where T is the duration of the adiabatic passage |5l| . In 
addition we note that in the adiabatic limit, the excited 
state |c) does not get populated during the entire process, 
so that spontaneous emission plays no role. 

In the second step of our unitary restoration of the 
quantum information to the original basis states we co- 
herently split the electronic population associated with 
the motional state |4, 0), in the superposition of Eq. (57), 
without affecting the state |2, 2). This is accomplished 
through a unitary time evolution governed by the Hamil- 
tonian 

H split = hg(t) {&l a x - a\ a y ) <g> \a)(b\ + H.c. , (59) 



whose generation through two pairs of lasers that drive 
a stimulated Raman transition between the electronic 
states | a) and follows closely the manner in which 
the Hamiltonian given in Eq. (|32| ) was constructed. The 
first pair of Raman lasers is arranged so that the wavevec- 
tor difference ^fc'- 1 -' is aligned with the x axis, whereas, 
for the second pair of Raman lasers, we assume that the 
wavevector difference Sk^ is now aligned with the y axis. 
Following Eq. @, this gives rise to the Hamiltonians 

H^' 2 ) = -HgWtyexp [-??72] |o)(6| 

<g){l-77 2 al a x , y } +H.c. , (60) 



for the first and the second pair of Raman lasers, re- 
spectively. Here, we have further assumed the Lamb- 
Dicke limit ]27|^^| , and rjx — rj y = 77, for the Lamb- 
Dicke parameters, which can be accomplished through 
an appropriate choice of |(5fc| W an d |<5fc|' 2 ' as seen from 
Eq. (pll). Given that the laser phases are arranged such 
that !P^(t) = —g( 2 \t), the dynamics generated by the 
combination of the two pairs of Raman lasers is described 
by the Hamiltonian in Eq. (|59|), and the coupling con- 
stant is given by Eq. ([sBJ) . The resulting time evolu- 
tion is identical to the one in Eq. (|36|), but where now 
X = \g\{d x d x ~ a y dy). With the particular choice 
A = 7r/16|<?|, for the generalized pulse area in Eq. (|2S|), 
the Raman-induced dynamics leave the ion in the state 



|^ = i{(c + + O|4,0)®(|a> + |6>) 
+ (c+-c_)\/2|2,2)®|&)}, 



(61) 



after the second step in our unitary restoration. 

In the third and final step, we complete the transfor- 
mation in Eq. ( |53|) through a second adiabatic transfer 
process which is similar to the one that we have employed 
in the first step. Starting from the state given in Eq. (|6l|), 
we recombine the electronic population in the state \b) by 
adiabatically transferring the population from the elec- 
tronic state I a) to the state \b), such that the motional ex- 
citation is simultaneously transferred from the x into the 
y direction. This is done through an adiabatic passage in 
counter-intuitive pulse sequence as described above but 
where we now assume that the laser which generates the 
pump pulse is aligned with the x axis and tuned to the 
fourth red sideband of the \a) |c) transition. The laser 
that generates the Stokes pulse is tuned to the fourth red 
sideband of the \b) ^ |c) transition and aligned with the 
y axis. As discussed in Appendix [B|, this leads to the 
Hamiltonian 



,4 

Hcomb = -hg a (t) al ®|a)(c| 
-hg b (t) a)y® \b)(c\ 



H.c. 



(62) 



which has the important property that U com b = 
e -iH comb t/n = 1 on the states |4,0) <g> \b) and |2,2) <g> |6>. 
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The only component in the superposition state of Eq. ( |6l| ) 
which is affected by the adiabatic transfer is the state 
|4, 0) (g> \a). This state coincides with a component of the 
dark state associated with the above Hamiltonian 



IVWny (*)> dark 

a(t)\n x + 4, n y 



® | a) - ( 9(t)|n !C ,nj / + 4) <g> |6> 



(63) 



for n x 



0, before the pump pulse is turned on. 



As before, the time dependence of the dark state is given 
through the quantities a(t) = gl(t)[(n y + 4)!/n y !] 1 / 2 and 
(3(t) = g* a (t)[{n x + Ay./nJ} 1 / 2 , so that after the Stokes 
pulse is turned off, the state |4, 0) ® \a) has been adiabat- 
ically transferred to the state |0, 4) (g> \b). The adiabatic 
process generated by the pump and the Stokes pulse then 
leaves the ion in the product state \tp) = ®\b), 
where the motional state 



+ (c+ ~ 



r-c_)(|4,0> + |0,4)) 
c_)V2|2,2)}, 



(64) 



and where we have chosen the two lasers beams to be 
7r out of phase so as to cancel the accumulated phase 
acquired by |0, 4) ® \b). The final state in Eq. ( |6~| ) is iden- 
tical to the initial motional state in Eq. (||) , and thus the 
stable motional basis states given in Eq. ( p2| ) have been 
completely restored. 

To summarize the various steps in the above proposed 
restoration of the qubit after the detection of a quantum 
jump associated with the C x channel: we first adiabati- 
cally increased the excitation number of the ion's motion 
along the x direction using two time-delayed laser pulses. 
Secondly we coherently split the electronic population 
associated with the motional state |4, 0) into the super- 
position | a) + | b) using a unitary transformation which 
we generated through a stimulated Raman process in- 
duced by two pairs of laser beams. The third step of the 
restoration sequence then served to establish the super- 
position |4, 0) + |0, 4) in the motional basis states and at 
the same time disentangled the ion's motional and elec- 
tronic degrees of freedom. This final step was achieved 
through a second adiabatic transfer process generated by 
two time-delayed laser pulses. 

To complete our analysis of the time scales associated 
with our stabilisation scheme, we examine here the adia- 
batic transfer processes involved in the unitary inversion 
of the detected quant um ju mps. From our discussion 
at the end of Section IV A we can estimate the dura- 



tion of the laser pulse that generates the intermediate, 
resonant step which led us from Eq. (j57]) to Eq. (^lj) , to 
be Tl ~ 8fis. The minimal duration of the adiabatic 
transfer steps is limited by the adiabaticity condition 
in Eq. (58). This is essentially controlled through the 
magnitude of the coupling constants \g a \ and \g \. We 
consider here the second adiabatic transfer process since 



this involves the excitation of the ion's fourth vibrational 
sidebands which is highly suppressed through the small- 
ness of the Lamb-Dicke parameters T) Xjy . More specifi- 
cally \g a ,b\ = vi, y \9a,b\/^-i as seen from Eq. (B6) where 
\9a\ (\db\) gives the resonant coupling strength for the 
dipole transition between the states |o) (\b)) and |c). The 
magnitude of the resonant coupling strengths is in turn 
limited by the off-resonant excitation of other sidebands 
and in particular the carrier transition which we have 
neglected in making the vibrational rotating wave ap- 
proximation. The mathematical condition for this limit 
is \g a ,b\/4vx,v < 1 fH], and for v XtV /2-K = 10MHz, im- 
poses \g a ,b\/^ <C 40MHz, on the magnitude of the reso- 
nant coupling strengths. We assume \g a ,b\/^ = 15MHz, 
which for r\ x , y = 0.2, leads to \g a ,b\ = 1kHz, for the 
coupling strength corresponding to the excitation of the 
fourth vibrational sidebands. With this result the adia- 
baticity condition in Eq. H imposes the limit 



T L > 320^s , 



(65) 



on the duration of the pulses that generate the adiabatic 
transfer and where we have assumed the pulse shape 
given in Eq. (|4l]) . The duration of this process may 
be shortened for larger values of the resonant coupling 
strengths |<7 a ,&|i which would however require a simul- 
taneous increase in the trap frequencies so as to avoid 
the off-resonant excitation of other than the chosen side- 
bands. 

Following our assessment for the duration of the 
various processes involved in our stabilisation scheme 
(Eqs. (ft2]), ( |5l| ) and (|65"|)), we can estimate the total du- 
ration of steps (i) and~(ii) indicated in Fig. [j] to be of the 
order of 1ms, with currently available technology. With 
this result we can give a figure of merit for the perfor- 
mance of our scheme. We assume r = 10ms, for the 
period with which we interrogate the dissipative system 
evolution and 7 = 0.1s -1 for the decoherence rate. With 
these values the probability for two quantum jumps to 
occur during the time interval r equals the probability 
for a single quantum jump to occur during the interroga- 
tion process and is of the order of 10~ 2 . Our stabilisation 
scheme then suppresses the rate of decoherence by two 
orders of magnitude. We note that this is not a funda- 
mental limit and can be improved through reducing the 
duration of the operations that constitute step (i) and 
(ii) in Fig.|l|. As discussed above this would primarily 
require an increase in the trap frequencies v x . y , beyond 
current laboratory values. 



VI. CONCLUSION 

To conclude, we have considered the active stabilisa- 
tion of a qubit which is stored in the bosonic motional de- 
grees of freedom of an ultra-cold ion and which is suffer- 
ing a particular type of motional decoherence. Through 
a quantum-trajectory analysis of the dissipative system 
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evolution associated with a finite-temperature master 
equation we were able to uncover particular basis states 
with which to encode the qubit and which obey the nec- 
essary criteria for the effect of a single quantum jump and 
the absence of quantum jumps associated with this form 
of decoherence to be unitarily reversible. We found that 
there exists a duality between the different unravelings 
of the master equation and the form of the stable basis 
states associated with each unraveling. Based on this du- 
ality we chose a form for the basis states such that the 
signature of the quantum jumps can be extracted through 
two consecutive binary interrogations. These interroga- 
tions determine if a single jump has taken place and if 
so, which type of jump. The final stage of the interro- 
gation involved the reading out of the information held 
in the electronic states. We proposed a method which, 
through coupling the ion to a low-Q optical resonator, 
can transfer this electronic information to the environ- 
ment with very little disturbance to the motion. Finally, 
the unitary inversion of the detected quantum jumps re- 
quired two adiabatic processes and one stimulated pro- 



cess. The entanglement (Section IV A), and unitary in- 
version (Section [v|), steps utilized, at most, two pairs 
of Raman lasers. This should be experimentally feasi- 
ble. The repair operation described here is completely 
unitary, and does not make recourse to any non-unitary 
probabilistic operations. Thus it should be ideally suited 
to operate on single systems such as the ions in a quan- 
tum computer. In addition to being stable against the ef- 
fects of a zero-temperature bath, we have shown that the 
particular encoding found is also stable against the effects 
of a thermal bath. This result might prove very signifi- 
cant as thermal noise plays a large role in the systematic 
sources of decoherence in trapped ions [fl3|| . However, 
we have restricted our discussion regarding the quantum 
jump detection and inversion processes to the case of a 
zero-temperature bath. The case of a thermal bath will 
be better formulated when the separate issue of the gen- 
eration of any given motional unitary operator has been 
solved. The entire process of restoration was made pos- 
sible by the a priori knowledge of the type of dissipation 
present. In a typical experiment, there will be many 
varied sources of noise and for the above stabilisation 
to be useful, information concerning the various sources 
decoherence must be experimentally obtained. Finally, 
it would be very interesting to learn whether the above 
formal analysis concerning the existence of stable bases 
extends to more complicated types of dissipation. The 
practical implementation of such extensions will, how- 
ever, pose a formidable task. 
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APPENDIX A: PULSED STIMULATED RAMAN 
TRANSITIONS 

The transition between the two ground state hyperfine 
levels | a) and \b) can be driven by two lasers connecting 
these two states to a common excited state |c), in a stim- 
ulated Raman scheme |}3|,[34]], as indicated in Fig.|2|. We 
briefly review this here for the case of a pulsed transition 
and derive som e of th e notation necessary for our discus- 
The stimulated Raman transition 



sion in Section IVA 



is driven by the electric field 

E(t)=E a (t)e- i &~*- 
+ E b (t)e-^ 



"*] + 



(Al) 



where E a {t) — E a f{t), and 23&(t) = Ebf(t) are the elec- 
tric field amplitudes of the two lasers and f(t) is a dimcn- 
sionless pulse shap e, wh ich we assume to be the same for 
both lasers. In Eq. ( Al) we have denoted the frequency of 
the laser driving the |z) |c) transition by uj C i, and the 
corresponding wavevector is given by where i = a,b. 
To generate the stimulated Raman transition between 
the two states \a) and |6) without populating the upper 
level |c), we assume both lasers to be far detuned from 
the excited state | c) and denote the corresponding detun- 
ing by Ai = (uj c —LUi) — uj C i. We are considering a resonant 
two-photon transition between the states \a) and \b) so 
that 

A a = A 6 = A , (A2) 

for the detuning of the two lasers from the excited state 
|c). In the dipole approximation, and after performing 
the optical rotating wave approximation, we can adiabat- 
ically eliminate the excited state |c) under the conditions 
|A| ^> \g a (t)\,\gb(t)\, as shown in Q. Here we have 
defined the dipole coupling constants gi{t) = f(t)gi — 
f(t) (i\p\c) Ei I H, and p is the dipole operator. After the 
adiabatic elimination we obtain the Hamiltonian 



H 



-hS a (t)\a)(a\-hd b (t)\b)(b\ 



-hg(t) exp 



E 



(r£+r$)/2]\a){b\ 
(-ir)j) m i +n * 



j—x,y 11 j ,mj =0 



rrijlrij 



-t ' 



exp [—i Vj (rij — rrij) t] + H.c. 



(A3) 



in the interaction picture of Hq, which was given in 
Eq. (g). Here we have introduced the Raman coupling 
constant 
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g{t) = f{t? 9 ; 



(A4) 



where g — g a g b /2A, and defined the Lamb-Dicke param- 
eters r\ x — Axo S k x and r\ y = A?/o Sk y , through the x and 
?/ component of the wavevector difference Sk = k a — k b , 
and the width of the motional ground state along the 
x and y axes as given by Axo = (h/2h> x m) 1 ^ 2 and 
Ay = (h/2v y m) 1 / 2 . 

The stimulated Raman transition has the disadvantage 
that it gives rise to optical Stark shifts which are given 
by the terms 



S a (t) 



A 



S b (t) = 



A 



(A5) 



and which are time dependent in the case of a pulsed 
transition JmJ. If these shifts are equal, i.e. for \g a (t)\ 2 — 
\gb(t)\ 2 , the relative energy shift between the states \a) 
and 1 6) is constant and zero, and the optical Stark shifts 
will merely contribute as an overall phase to t he dynam- 
ics generated by the Hamiltonian in Eq. (A3). We will 
assume 5 a (t) = 5b(t) and drop the contribution of these 
optical Stark shifts, thereby not explicitly keeping track 
of overall phases which are not important to us here. This 
situation can be experimentally realized by varying the 
relative intensity of the two Raman pulses ||. Alterna- 
tively, it has been shown that appropriately chirped laser 
pulses can compensate for these time varying Stark shifts 



In our discussion in Section IV A we are interested in 



the situation where the stimulated Raman transition de- 
scribed by the Hamiltonian in Eq. ( [A3] ) is only sensitive 
to the motion of the ion along one of the principal axes. 
This can be realized by arranging the two exciting laser 
beams so that their wavevector difference 8k = k a — k b , is 
aligned wi th that principal axis. The Hamiltonian given 
in Eq. (A3) then simplifies to 



H 



-%(i)exp [-7^/2] \a)(b\ 



' E 

n.m— 



(-%) 



m-fn 



)J 71. 1 



x exp [— % Uj (n — m) t] + H.c. 



(A6) 



where j = x for the case of aligning 6k_ with the x axis, 
and j — y for the case of aligning (5fc with the y axis. 

A further simplification arises in the low excitation 
regime |42[ ], and if the pulse turn-on and off is slow 
when compared to the trap frequencies v x and v y . We 
can then perform the vibrational rotating wave approx- 
imation p8[, and consider only the resonant terms in 



Eq. (A6), to obtain 



H = -Hg(t)exp[-rfi/2] \a){b\ 



>E 

n=0 



n\n 



H.c. 



(A7) 



where j = x,y as explained above. Without induc- 
ing transitions between the ion's vibrational levels, this 
Hamiltonian is sensitive to the ion's motional state 
through its dependence on the creation and annihilation 
operators for the ion's vibrational excitation. 



APPENDIX B: ADIABATIC PASSAGE 

To effect the unitary inversion of the decoherence pro- 
cesses associated with the decay channels C x and C y , we 
make use of adiabatic passage techniques which have pre- 
viously been studied in the context of coherent popula- 
tion transfer pi] , where atomic population is coherently 
transferred via Raman transitions induced by a pair of 
overlapping time-delayed laser pulses. These techniques 
rely on the existence of a dark state of the corresponding 
Hamiltonian. In Section ^ we employ adiabatic passage 
processes that are based on the dark states associated 
with the two types of Hamiltonians 



H, 



dark 



-hg a (t) a^ y a (g) \a)(c\ 
-hg b (t) a x K " ® \b)(c\ +H.c. 



and 



Hdark = -%o(*) al a ® \a)(c\ 
-hg b (t) al Kb ®\b)(c\ 



H.c. 



(Bl) 



(B2) 



which describe the excitation of the A system depicted 
in Fig. ^| with two laser pulses that drive specific vibra- 
tional sidebands of the dipole transitions \a) \c) and 
| b) <^> |c), respectively. Here, we briefly describe the de- 
tails of this excitation. Vogel and de Matos Filho [EEl , 
have shown that the excitation of the dipole transition 
\i) |c), (where i = a,b) with an electric field 



E i {t)=f i (t)E i e- i ^- r --^ + 



c.c. 



(B3) 



which is tuned to the ftth red sideband of the transition 
is described by a nonlinear K-quantum Jaynes-Cummings 
model in the form 



H 



-fift(*)exp[-^/2] \i)(c\ 



1 E T) \( n j Ki i "./ 

n=0 v ' 



H.c. 



(B4) 



where j = x, for the case of aligning the wavevector fc 4 
with the x axis, and j — y, for the case of aligning with 
the y axis . Here we have defined the dipole coupling 
constants <?;(£) = fi{t)gi = /,(*) {i\p\c)Ei/h, where p is 
the dipole operator and the dimensionless quantity fi(t) 
gives the time dependence of the electric field amplitude 
as in Eq. ( |B3| ) . The Lamb-Dicke parameters are given 
by f]x = Ax^jfcJ and r\ y — Ay |fcil f° r the case of align- 
ing the wavevector with the x and y axis, respectively. 



In the Lamb-Dicke limit r] x ^ y <§; 1 



and we retain 
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only the leading order term in the Lamb-Dicke parameter 
to obtain 



H = -h~ 9l (t)e 



-ti?/a (-^i)" 



a] ®|i)(c|+H.c (B5) 



With this result, the pulsed excitation of the A system 
with two las er b eams is described by the Hamiltonian 
given in Eq. (Bl) if the laser that drives the \a) \c) 
transition is aligned with the y axis and tuned to the 
K a th red sideband, and the laser that drives the \b) |c) 
transition is aligned with the x axis and tu ned to the 
Kbth red sideband. The Hamiltonian in Eq. ( |B2[ ) is ob- 
tained for the same sideband detunings but where the 
laser that drives the \a) \c) transition is aligned with 
the x axis, and the laser that drives the \b) O |c) tran- 
sition is aligned with the y axis. The coupling constants 
in Eqs. (Bl) and (B2) are given by 



g i {t) = h{t)g i = g i {t)^[-ri 1 i /2 



{-ir]i) K * 



(B6) 



where i — a, b, and the Lamb-Dicke par am eters r\ a ,b — 
t]y tX in Eq. (Bl), and 7] a j, — r) x , y in Eq. (B2). The time 
delayed pulses are characterized by f a (t) and fb(t) which 
give the dimensionless pulse shape. 
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